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Abstract. - The non-equilibrium structural and dynamical properties of flexible polymers confined 
in a square microchannel and exposed to a Poiseuille flow are investigated by mesoscale simulations. 
The chain length and the flow strength are systematically varied. Two transport regimes are 
identifled, corresponding to weak and strong confinement. For strong confinement, the transport 
properties are independent of polymer length. The analysis of the long-time tumbling dynamics of 
short polymers yields non-periodic motion with a sublinear dependence on the flow strength. We 
find distinct differences for conformational as well as dynamical properties from results obtained 
for simple shear fiow. 
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Introduction. — Confinement fundamentally alters 
the properties of dilute polymer solutions - compared to 
bulk behavior - when either the polymer radius of gyra- 
tion is on the order of the characteristic dimensions of its 
proximity [1-4], and/or an external field is applied, e.g., 
a shear or pressure field. In the first case, geometrical 
constraints lead to a stretching of the polymer parallel to 
the surfaces [1-6] . In the second case, in addition to flow- 
induced deformations, polymer-surface hydrodynamic in- 
teractions determine the polymer dynamics and leads to, 
e.g., cross-stream migration [7-12]. The migration effect 
has been studied intensively for DNA-like molecules by 
computer simulations and the relevance of hydrodynamic 
interactions has been confirmed. Experimental studies of 
individual DNA molecules in steady shear flow by fluores- 
cence microscopy have provided a wealth of information on 
single polymer dynamics [13-16]. In particular, large con- 
formational changes have been revealed due to tumbling 
motion [15,17,18]. Similar studies for flexible polymers 
confined between surfaces or in narrow channels have not 
been conducted so far, however, a similar complex dynam- 
ics can be expected. Understanding of single polymer be- 
havior is of paramount importance for the emerging tech- 
nology of microfluidic devices. Insight into the detailed mi- 
croscopic conformational, dynamical, and transport prop- 
erties of polymers, e.g., DNA, will help in the conception 
and design of such devices. Moreover, such studies will 
contribute to the understanding of the transport proper- 
ties of biological macromolecules through blood vessels. 



The proper account of hydrodynamic interactions is es- 
sential in simulation studies of fluid flows in channels as 
is emphasized by the appearance of cross-streamline mi- 
gration. Recently developed mesoscale simulation tech- 
niques, such as Lattice Boltzmann simulations [12, 19], 
Brownian dynamics simulations with a hydrodynamic ten- 
sor [4], and multi-particle-collision dynamics [20,21] (also 
called stochastic rotation dynamics), are well suited for 
simulations of microchannels flows and are able to bridge 
the length- and time-scale gap among the solvent and 
solute degrees of freedom thereby taking hydrodynamic 
interactions adequately into account. In this letter, we 
will present results for the conformational and dynamical 
properties of polymers confined in a square channel and 
exposed to a Poiseuille flow by mesoscale computer simu- 
lations. Both, the polymer length as well as the pressure 
gradient are systematically varied. The considered poly- 
mers fall into the crossover regime from weak to strong 
confinement. Experiments [11] and simulations [8] predict 
a distinct different flow behavior in the two limits, which 
is conflrmed by our simulations. The strong confinement 
regime has received less attention so far, it is however most 
relevant if microchannels arc intended to be used for DNA 
characterization. In addition, we discuss chain orientation 
and tumbling dynamics, aspects which have not been ad- 
dressed in previous studies. 

Simulation method, model. We use a hybrid 
simulation approach to study the properties of flexible 
polymers in flow, where molecular dynamics simulations 
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Fig. 1: Conformations of a polymer of length A'^ = 160 in a 
channel of width L — 15a for weak (Peclet number Pe^ — 10, 
top two images) and strong flow [Pe^ = 100, bottom). The 
latter yields close to fully stretched polymers. 

(MD) are combined with the multi-particle-collision dy- 
namics (MPC) method for the solvent [22,23]. MPC is 
particle based and proceeds in two steps. In the stream- 
ing step, the solvent particles of mass m move ballistically 
for a time h. In the collision step, particles are sorted 
into the cells of a cubic lattice of lattice constant a and 
their relative velocities, with respect to the center-of-mass 
velocity of each cell, are rotated around a random axis 
by an angle a. For every cell, mass, momentum, and 
energy are conserved. The algorithm is described in de- 
tail in refs. [22,24-27]. The fluid is confined in a square 
channel of side length L and periodic boundary conditions 
are applied along the channel axis (i||). No-slip bound- 
ary conditions are imposed by the bounce-back rule, as 
described in Ref. [28], and flow is induced by a gravi- 
tational force {F = mg) acting on every fluid particle. 
The linear polymer is comprised of mass points of mass 
M, which are connected by a strong harmonic potential 
with an equilibrium bond length b. Excluded-volume in- 
teractions are taken into account by the well-known 12-6 
truncated Lennard- Jones potential with the parameters a 
and e [29]. Since we consider pressure-driven flows, no 
gravitational force acts on the polymer. The interaction 
of a polymer with the solvent is realized by inclusion of 
its monomers in the MPC collision step [30]. Extensive 
studies of polymer dynamics confirm the validity of this 
procedure [20,26,29,30]. Specifically, we employ the pa- 
rameters a — 130°, h — O.lr, with r = a/ ma?- /ksT (fcs is 
Boltzmann's constant and T the temperature), the num- 
ber of particles per cell p = 10, M = mp, b = a ~ a, 
the fluid mass density g — pm/a^, ksT/e = 1, and the 
time step in MD simulations /imd = 5 x 10^'^r. Poly- 
mers with the monomer numbers N = 20, 40, 80, and 160 
are considered. To minimize self-interactions of periodic 
images, the corresponding channel lengths are L|| = 40a, 
60a, 100a, and 170a, respectively. Simulations of the pure 
solvent confirm that the generated velocity profiles agree 
with the solution of the Navier-Stokes equations for the 
considered geometry. 

Results. — We expect a pronounced different behav- 
ior of polymers in Poiseuille flow for radii of gyration ei- 
ther smaller or larger than the width of the channel [11]. 
Scaling considerations based on a blob model [1] and sim- 
ulations [4] yield an equilibrium stretch of a flexible poly- 
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Fig. 2: Radial monomer distribution functions for the Peclet 
numbers Pe^ = 0, 10, 50, 100, 200 (top to bottom at r/rgo = 
1.5), A'^ = 40, and L/a = 15. Inset: widths of the distribution 
functions for the chain lengths TV = 20 (squares), A'' = 40 
(diamonds), TV = 80 (bullets), and TV = 160 (triangles). 

mer in a microchannel when its bulk radius of gyration 
Vgo ^ bN^^^ exceeds the lateral channel dimension, i.e., 
for ^ L. Here, the channel width becomes the char- 
acteristic length scale and the blob relaxation time, which 
is determined by L, will be the relevant time scale [4,8]. 
For <SC L, the conformations are unaffected by con- 
finement and the characteristic length and time scales are 
given by the bulk values of the polymer. The ratios 2rgo / L 
for the considered polymer lengths and the channel width 
L — 15a cover the range 0.36 — 1.3, corresponding to the 
crossover regime between weak and strong confinement. 
The snapshots of fig. 1 illustrate the preferred average 
alignment of a polymer along the flow direction. 

Figure 2 provides an example of the radial monomer 
distribution for the polymer length TV = 40 and various 
Peclet numbers. Here, we defined the Peclet number Pe^ 
as Pe^ = ITl, with an effective shear rate 7 — gLg/{2ri) 
and the relaxation time tl — ri{L/2)^ /{3kBT) {j] is the 
solvent viscosity) [31], where the blob diameter has been 
replaced by the channel width. A similar approach has 
been adopted in ref. [8]. The distribution functions are 
normalized such that Pm{r / r gQ)r / r'^f^dr — 1. With in- 
creasing Peclet number, the distribution Pm increases at 
r/rgo ~ 1 due to the stretching of the polymer by the flow. 
For even large Pe^ , migration sets in and the maximum 
of P„i shifts to smaller radii, as already observed in num- 
ber of other studies [7-10,12]. The widths of the radial 
distribution functions 

/"OO 

= / r^P^{r)dr (1) 

JQ 

for the various chain lengths are shown as inset in fig. 2. 
For small Peclet numbers, the widths are independent of 
Pe^ and decay at large Pe^. There seems to be a weak 
dependence on chain length only. The widths start to de- 
crease at approximately the same Peclet number. This 
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Fig. 3; Radial monomer distribution functions for the channel 
widths L/a = 15, 20, 30, and 45 (left to right). The polymer 
length is A'^ = 40 and gr^/a — 0.01, i.e., the curvature of 
the flow profile is independent of L. Inset: position of the 
maximum as function of the channel width. The straight line 
represents the power law (L/rgo)*^^ . 

is consistent with our assumption that the blobs deter- 
mine the polymer properties. The approach to a uni- 
versal limiting curve with increasing polymer length is 
explained as follows. In general, P,„ exhibits the de- 
pendence Pm{b, r, L, N, g), and correspondingly the width 
= (r2) = W^{b,L,N,g). Scaling aU lengths by L and 
using the above definition of the Peclet number yields the 
dependence W'^/L'^ = f{rgo/L,bN/L,Pe^). For strongly 
confined polymers, the width is determined by blobs - and 
thus by L - rather than by the polymer length, hence, 
we find W'^/L'^ — f{Pe^), i.e., it depends on the Peclet 
number only. According to current understanding, mi- 
gration is caused by polymer-surface hydrodynamic inter- 
actions [7-10, 12]. Since hydrodynamic interactions are 
long-range in nature, it is difficult to unravel the system 
size dependence of migration. To study the influence of 
the channel width on migration, we performed simulations 
for various channel widths L at the same curvature of the 
flow proflle {gr'^/a — 0.01) and polymer length N = 40. 
As shown in fig. 3, cross-streamline migration is well pro- 
nounced for all L, with a shift of the maximum of the 
monomer distribution to larger radii with increasing L. 
The position of the maximum (r^) exhibits the chain- 
length dependence ~ L^^^ as shown in the inset of 
fig. 3, i.e., migration is more pronounced for less strongly 
confined polymers. The changes in polymer conforma- 
tions are refiected in their radii of gyration, which are 
displayed in fig. 4. Evidently, short polymers {N < 40) 
experience only minor conformational restrictions at equi- 
librium, whereas longer polymers are stretched along the 
channel and squeezed in the transverse direction. Beyond 
a certain Peclet number Pe^, flow induces or enhances 
stretching parallel to the channel and transverse shrink- 
age. This Pe^ value depends only very weakly on chain 
length. 
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Fig. 4: Mean square radii of gyration parallel (squares) and 
transverse (bullets) to the channel for the chain lengths A'' — 
20, 40, 80, 160 (bottom to top, squares; top to bottom, bullets) 
and L — 15a as a function of the Peclet number Pe^ . 

Figure 5 shows the component of the mean square 
radius of gyration along the channel axis for various 
chain lengths. Here, a different kind of Peclet num- 
ber is employed, which is deflned as Pe = jrc, with 
To = rjrgg/ {3kBT) [31], i.e., the polymer relaxation time 
is used instead of a blob-related time. The Peclet num- 
ber is then chain- length dependent. The shorter polymers 
seem to exhibit universal behavior for Peclet numbers up 
to Pe « 5. For larger Pe, flnite polymer-size effects play 
a role. However, the longer polymers show a different be- 
havior. This confirms our picture of the various involved 
length and time scales. The polymers of lengths < 40 
are only weakly confined in the channel. The flow behav- 
ior is determined by the equilibrium radius of gyration and 
the longest polymer relaxation time. Instead of the effec- 
tive shear rate 7 = gLg/{2ri), we could use a value deflned 
in a different way, e.g., the local shear rate averaged over 
the channel width. For longer polymers {N > 40), the 
channel width deflnes the relevant length and time scales 
as long as the confinement is felt by the polymer. At high 
fiow velocities, the polymer migrates to the channel center 
and is no longer in contact with the surfaces (cf. fig. 2). 

So far, we have discussed conformational properties 
which are averaged over the microchanncl. Due to the 
parabolic velocity profile, however, the local shear rate 
changes linearly with the radial distance. Thus, we can de- 
fine a local Peclet number by Pe'^™ — ggr^g^rcm/ {iksT), 
where rem is the radial chain center-of-mass position. Fig- 
ure 6 shows the dependence of the mean square radius of 
gyration parallel to the channel as function of Pe'^" for 
the chain length = 40 and L = 15a. The various data 
sets cover the range of Peclet numbers Pe^ = 1 — 200. For 
small Pe"^™, i.e., in the central part of the channel we find a 
stretching of the polymer by the flow field which increases 
with increasing field strength Pe^. Thus, in the central 
part, the channel flow is different from a simple shear flow 
- the extended polymer never experience a linear veloc- 
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Fig. 5: Mean square radii of gyration parallel to the channel 
for the chain lengths A'^ = 20, 40, and 80 (bottom to top) and 
L — 15a as a function of the Peclet number Pe. 

ity profile only. At larger radial distances, the polymer 
is further stretched, particulary for larger Pe^ [9]. More 
importantly, the various data sets seem to approach a lim- 
iting curve for sufficiently large rem- However, this curve 
deviates from the theoretical dependence valid for linear 
shear flows, at least for the length range covered by our 
simulations. The transverse radius of gyration exhibits 
a similar deviation from the theoretical prediction. This 
leads us to the conclusion that the non-equilibrium struc- 
tural properties of flexible polymers in a Poiseuille flow 
are even locally different from those of a polymer in lin- 
ear shear flow at least as long as the characteristic chain 
dimensions are on the order of the channel width. The 
reason seems to be the variation of the relevant relaxation 
times as Pe^ and the center-of-mass position change. 

In addition, we analyze the orientational distribution 
function of the end-to-end vector. Since the system is 
almost rotational symmetric with respect to the chan- 
nel axis, only the distribution function Pa{d), where 9 
(0 < < tt) is the angle between the end-to-end vector and 
the flow direction, needs to be studied. As shown in fig. 7, 
the distribution is almost independent of the angle 9 at 
small Peclet numbers. The deviation from the bulk value 
1/2 is caused by the confinement of the polymer. (The 
distribution is normalized such that Pa{9) sin 9d9 = 1.) 
As expected, there is no preferred angle; the average is 
(9) — 0. The shape of the distribution strongly depends on 
the Peclet number. For large Peclet numbers. Pa decays 
algebraically as 9~'^. A similar algebraic decay is found 
for the angle between the end-to-end vector and its pro- 
jection onto the shear plane of a flexible polymer in shear 
flow [16, 18]. However, the distribution Pa also reflects 
the fundamental difference of polymers in Poiseuille and 
shear flow. For shear flow, there exists a regime of shear 
rates, where the distribution is Gaussian, only for large 
shear rates follows the algebraic decay [16,18]. In fig. 7 no 
Gaussian regime is present. Because of the nonlinear flow 
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Fig. 6: Mean square radii of gyration parallel to the channel 
for the Peclet numbers Pe^ = 1, 5, 10, 20, 30, 50, 100, 200 
(left to right) as a function of the Peclet number Pe'^'", where 
Pe*^'" ~ rem- The solid line is an analytical result for a chain 
of the same persistence length in simple shear flow [18]. 



profile, an algebraic decay is not necessarily surprising, 
compared to simple shear, where this fact is explained by 
higher order correlations and the importance of the whole 
time history. Nevertheless, the similarity is striking and 
might be traced back to the same origin. The radial width 
of the monomer and center-of-mass distributions decrease 
with increasing flow strength (cf. fig. 3) and a polymer 
spends more time in the region of higher axial fiuid veloc- 
ity [8]. This results in an increased polymer transport ve- 
locity compared to the average solvent velocity, as shown 
in fig. 8. For a wide channel r^o <C L, the ratio of Vcm/vs 
is close to unity at equilibrium. Conformational restric- 
tions for <g; L lead to a narrowing of the radial distri- 
bution functions [11] accompanied by increased polymer 
flow velocities even at very weak imposed flows. This ex- 
pected transport behavior is shown in fig. 8. At small 
Peclet numbers Pe^ < 50, we find - except for = 20 
- center-of-mass velocities which depend only weakly on 
Pe^, caused by confinement. Only for the polymer of 
length A = 20 reduces the velocity ratio to unity at equi- 
librium. With increasing Peclet number, the average poly- 
mer velocity increases, starting at about the same Peclet 
number Pe^ for all chain lengths, and the ratio reaches 
almost its maximum value of two. We obtain a chain- 
length dependence of the velocity ratio for short polymers 
only, which is rather a consequence of the weaker confine- 
ment of the shorter polymer than a true molecular weight 
dependence. The ratios for the longer polymers seem to 
approach a universal limiting curve, which is consistent 
with the assumption that local properties on a blob scale 
rather than global polymer properties determine their be- 
havior. As a consequence, no separation of highly confined 
DNA according to molecular weight is possible. Polymers 
in simple shear flow exhibit large conformational changes 
due to tumbling motion, i.e., they stretch and recoil in the 
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Fig. 7: Probability distribution of the angle between the end- 
to-end vector and the flow direction for the Peclet numbers 
Pe^ = 1, 5, 10, 30, 50, 100 (bottom to top) and the channel 
width L — 15a. Only half of a period is shown. The solid line 
represents the power-law dependence Pa ~ . 




Fig. 8: Ratio of the polymer center-of-mass velocity Vcm and 
the average solvent velocity Vs for the chain lengths A*' = 20 
(blue), 40 (red), 80 (green), and 160 (orange) and the channel 
width L = 15a. 



com'se of time [15-18]. Visual inspection of the polymer 
dynamics in Poiseuille flow reveals similar large conforma- 
tional fluctuations (cf. fig. 1). Tumbling is characterized 
by a typical time denoted as tumbling time tx- There 
are various ways to extract such a characteristic time, al- 
though the dynamics is non-periodic. We tried various 
strategies and found the following approach most useful. 
We first determined the distribution of times between con- 
secutive crossings of the end-to-end vector with the plane 
perpendicular to the flow direction, i.e., 9 = ±7r/2. As 
shown in the inset of fig. 9, this yields exponentially de- 
caying distributions {Pt ~ e"*/*"^) at large times. The 
extracted tumbling times It are presented in fig. 9. A 
power-law fit yields the dependence tx ~ 

Sim- 
ulations and analytical calculations [15-18] yield the de- 
pendence tT ~ for flexible polymers in linear shear 
flow. Hence, the tumbling times of polymers in Poiseuille 
flow exhibit a weaker dependence on the strength of the 
external fleld than in shear flow. 

Conclusions. — We have analyzed the flow behavior 
of flexible polymers confined in a square channel with a 
side length comparable to the polymer radius of gyration. 
We have identified two distinct transport regimes were 
polymers behave similarly. For r^o ^ L the polymers are 
unaffected by the channel geometry and the characteristic 
length and time scales are given by the radius of gyration 
and the longest polymer relaxation time, respectively. In 
the strongly confined regime r^o ^ L, the polymer struc- 
ture and dynamics are determined by blobs with a size 
of the channel cross section. Polymers with r^p ~ L/2 
show a complicated and non-universal behavior. We ex- 
pect that the properties of polymers in a slit geometry are 
less severely affected by confinement than those of poly- 
mers in a channel. In the slit geometry, the polymer is 
always able to relax in the direction parallel to the sur- 
faces and hence the longest relaxation time is always re- 



lated to the polymer length. The analysis of the radial 
polymer deformation, the distribution of the orientation 
angle, and the tumbling times reveals differences among 
Poiseuille flow and simple shear flow for the strongly con- 
flned polymers. 
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